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We show that the QCD Dirac spectrum at finite chemical potential using a matrix model in the 
spontaneously broken phase, is amenable to a generic 2-dimensional effective action. The eigenvalues 
form a droplet with strong screening and plasmon oscillation. The droplet is threaded by a magnetic 
vortex which is at the origin of a Berry phase. For quarks in the complex or Dirac representation, 
the anomalons transport in the droplet of eigenvalues bear some similarities with that in droplets 
of composite fermions at half filling snggesting that the latters maybe Dirac fermions. 
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I. INTRODUCTION 

QCD breaks spontaneously chiral symmetry with a 
wealth of evidence in hadronic processes at low ener¬ 
gies m- First principle lattice simulations strongly sup¬ 
port that [2] . The spontaneous breaking is characterized 
by a large accumulation of eigenvalues of the Dirac oper¬ 
ator near zero-virtuality [3]. The zero virtuality regime 
is ergodic, and its neighborhood is diffusive [4]. 

The ergodic regime of the QCD Dirac spectrum is 
amenable to a chiral random matrix model In short, 
the model simplifies the Dirac spectrum to its zero-mode- 
zone (ZMZ). The Dirac matrix is composed of hopping 
between N-zero modes and N-anti-zero modes because of 
chirality, which are Gaussian sampled by the maximum 
entropy principle. The model was initially suggested as a 
null dynamical limit of the random instanton model [6]. 

QCD at finite chemical potential n is subtle on the 
lattice due to the sign problem [7]. A number of effec¬ 
tive models have been proposed to describe the effects 
of matter in QCD with light quarks [1]. Chiral random 
matrix models offer a simple construct that retains some 
essentials of chiral symmetry both in vacuum and mat¬ 
ter. For instance, in the chiral 1-matrix model finite /r 
is captured by a constant deformation of Gaussian ma¬ 
trix ensembles [siin]. In the chiral 2-matrix model the 
deformation with fi is also random [mill]. Chiral ma¬ 
trix models in matter were discussed by many nmg. 
Recently both a universal shock analysis [Hj and a hy- 
drodynamical description of the Dirac spectra were sug¬ 
gested m- 

In the ergodic regime the 1- and 2-matrix models ex¬ 
hibit the same microscopic universality for vanishingly 
small in the large volume limit m- This limit cor¬ 
responds to a weak non-hermitean deformation of the 
standard chiral matrix models and therefore preserves 
the underlying chiral symmetry of the coset space. It 
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follows the general strictures of the epsilon-expansion in 
chiral power counting |I6j . The microscopic universality 
for small fj? leads to new Leutwyler-Smilga sum rules m 
for the eigenvalues of the QCD Dirac operator that put 
first principle constraints on the effective Lagrangian ap¬ 
proach in matter as well QCD lattice simulations in the 
regime of small quark masses, small fj? and large vol¬ 
umes. 

In the first part of the paper, we will show that 
at finite fj- the distribution of Dirac eigenvalues in the 
complex plane maps onto a 2-dimensional Coulomb gas 
whose effective action is constrained by Coulomb' s 
law, the conformal and gravitational anomalies in 2- 
dimensions. These contributions are generic and go be¬ 
yond the specifics of matrix models at finite /r. The map¬ 
ping offers a physical framework through anomalies, for 
understanding first principle aspects of the QCD Dirac 
spectrum at finite /i with the hope of constraining fur¬ 
ther the effective Lagrangian approaches and current and 
future lattice simulations near the chiral point at finite 

and for large volumes. We will also show that the 
2-dimensional Coulomb gas of eigenvalues exhibits both 
screening and a plasmon branch. The latter translates 
to a diffusive mode in the stochastic evolution in chiral 
matrix models |I5j . The plasmon frequency sets an esti¬ 
mate for the relaxation time for the restoration/breaking 
of chiral symmetry at finite /r, solely through the stochas¬ 
tic re-organization of the low-lying modes of the Dirac 
spectrum. 

In the second part of the paper, we follow by noting 
that under adiabatic changes the complex eigenvalues at 
finite viewed as particles in the complex 2-plane, be¬ 
have as fermions for /3 = 2 due to the emergence of a 
Berry phase of tt. There is no Berry phase for /3 = 1 
or when the quarks are in the real representation. We 
use this observation to note that a dynamical droplet of 
quark eigenvalues for (3 = 2 and finite /r share some sim¬ 
ilarities with composite fermions at half filling suggest¬ 
ing that the latters are Dirac particles. We show that 
the anomalous charged transport contributions from an 
induced Wess-Zumino-Witten term are consistent with 
constituents of charge e = 1 and spin s = 1/2. We use 
these observations to derive novel effects for the compos- 
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ite fermions at half filling. 

This paper consists of the following new results: 1/ 
a generic anomalous effective action for the QCD Dirac 
eigenvalues as a droplet in the complex 2-plane of eigen¬ 
values with scalar curvature; 2/ a plasmon dispersion law 
in the bulk of the droplet with an improved estimate for 
the relaxation time for the breaking/restoration of chi¬ 
ral symmetry in QCD at finite fx; 3/ an identification 
of the eigenvalue droplet as a 2-dimensional Fermi liq¬ 
uid threaded by a magnetic vortex; 4/ a description of 
the anomalous transport on the Fermi surface through a 
Berry induced Wess-Zumino-Witten type term; 5/ a sug¬ 
gestion through geometry that composite fermions in the 
fractional quantum Hall effects at half filling are Dirac 
particles; 6/ two novel anomalous transport effects in 
composite fermions at half filling, caused by a rotation 
or a temperature gradient. 


II. 2-MATRIX MODEL 



FIG. 1: Eigenvalue distribution from a 2-matrix model. 
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The low lying eigenmodes of the QCD Dirac opera¬ 
tor capture some aspects of the spontaneous breaking 
of chiral symmetry both in vacuum and in matter. Re¬ 
markably, their fluctuations follow by approximating the 
entries in the Dirac operator by purely random matrix 
elements which are chiral (paired spectrum) and fixed by 
time-reversal symmetry (Dyson ensembles). At finite fx 
the Dirac spectrum on the lattice is complex [5D]. The 
matrix models at finite fx liiin] capture this aspect of the 
lattice spectra and the nature of the chiral phase tran¬ 
sition [niiaiia]. For a 2-matrix model, the partition 
function is [iniiii] 


Z[|i,^x,Nf,mf] = 

J g-a7VTr(BtB) 

xdetl 

A1 — ifxB^ —iruf 


( 1 ) 


for equal quark masses my in the complex representa¬ 
tion. Here A,B are valued, ix accounts for 

the difference between the number of zero modes and 
anti-zero modes. The vacuum Banks-Casher formula [3] 
fixes the dimensionfull parameter ^/a = Ig^gjo/n in terms 
of the massless quark condensate and the density of zero 
modes n = N/V 4 . Throughout, we will set the units 
using y/a —> 1. All canonical units are recovered by in¬ 
spection. 

The Dirac matrix in Q has ix unpaired zero modes 
and N paired eigenvalues iizj in the massless limit. The 
paired eigenvalues delocalize and are represented by 0. 
The unpaired zero-modes decouple. In terms of the fixed 
eigenvalues and large but finite N, Q reads unmi], 


The action is 

N N 

S[|3,^x;z] =-/3 \n\zf - (3) 

i<j=l i=l 


with 


W{z) =—a\a\z\ +—w{z) (4) 

and the quasi-harmonic potential r(;( 2 :) = \z\'^ —'^{z^+z‘^). 
Here a = /3(^ -I- 1) — 1 and ^ accounts for the difference 
between the number of zero modes and anti-zero modes 
in the ZMZ. We define 


A = ^ = l±if! = 1 

N/3 1 - t 2 2/r2 - p ^ 

with H = 1 / 1 % acting as a magnetic field with magnetic 
length Ib as we will suggest below. Throughout /3 = 2 
unless indicated otherwise. The 2-matrix model for /I = 
1,4 |21j is more subtle at finite p [la¬ 
in Fig. I^we display the distribution of eigenvalues fol¬ 
lowing from the 2-matrix model with A and B sampled 
from a Gaussian ensemble of 200 x 200 matrices with 
IX = Q and p = 0.3. The mean density in the droplet is 
po ~ vB/^'k with V = 1//3. The boundary curves follow 
from the analysis in |12j . The domain is an ellipse 



( 6 ) 


with semi-axes a^/2/2 = lir/lTT- The ellipse re¬ 
mains un-split with area A = 7ra+a_ = 27rZ^ for all values 








3 


of fjL. For the other quark representations with /3 = 1,4 
the joint distribution in the 2-matrix model is more sub¬ 
tle m- Throughout, will be assumed for /3 = 2, but 
most results extend to /3 = 1,2,4 for large N. 

For completeness we recall that the 1-matrix model 
corresponds to setting B = 1. In this case, the eigenvalue 
distribution forms a connected droplet in the z-plane for 
fj, < He and splits into two droplets symmetric about 
the real-axis for /r > fj-c, restoring chiral symmetry [3 
[9]. In the spontaneously broken phase, all droplets are 
connected and symmetric about the real-axis 

Throughout, ([^ is to be understood in the large N 
limit to allow for a course graining of the eigenvalue den¬ 
sity. We note that in this limit, the Bessel kernel in [10] 
is expanded. While the 2-matrix model is exactly solv¬ 
able in terms of the orthogonal polynomial method, our 
analysis of this model shows the emergence of a generic 
effective action for the complex eigenvalues of the Dirac 
spectrum at finite /r. Therefore our analysis encompasses 
the weak non-hermiticity range of the model. Besides chi¬ 
ral symmetry through the pairing of the complex eigen¬ 
values, the new guiding principles for the construct of 
this effective action are Coulomb law, the conformal and 
gravitational anomalies in 2-dimensions as we now detail. 


III. STATIC EFFECTIVE ACTION 

In this section we will re-write ([^ in terms of the ef¬ 
fective potential generated by the mean charged density, 
for a sufficiently dense ensemble of eigenvalues. This 
assumes that N is large to allow for a smoothening of 
the eigenvalue density over distances larger than the 
inter-level spacing but much shorter than the size of the 
eigenvalue droplet. While the 2-matrix model is solv¬ 
able [ini[TT], it is important to stress that its re-writing 
at large but fixed N, using an effective action unravels 
universal physical aspects of the eigenvalue droplet that 
are not restricted to the specifics of the model. We will be 
able to go beyond the strictures of chiral symmetry by 
making use of emergent anomalies when using the col¬ 
lective potential sourced by the smoothened eigenvalue 
density as we detail below. 

For the sake of generality, we will assume that the 
eigenvalue space is curved with a measure = gzzdzdz. 
The explicit form of the metric is not important for our 
general arguments. The curved eigenvalue space will al¬ 
low for the unraveling of two generic contributions to the 
effective action through the emergence of conformal and 
gravitational anomalies, that make the effective action 
construct more general than the 2-matrix model. Also, 
for bulk quantities a fixed eigenvalue curvature R is con¬ 
jugate to a fixed space curvature i ?4 « 1/i?. 

The re-writing of ( [^ w ill closely follow the effective ac¬ 
tion construction in ]22| for the Lauhglin states to which 
we refer for further details. In brief, on a curved eigen¬ 
value manifold of volume V 2 with a metric gzz and large 
N, the ensemble described by Q allows the change in the 


measure, from integrating over specific eigenvalues Zi to 
integrating over the eigenvalue density p{z). Specifically 

N 

i=l 

with the induced Boltzmann entropy in the exponent [181 
I22j . The mean density po be made explicit be¬ 
low. Since the spectrum is chirally symmetric, we reset 
{p{z) + p{—z))/2 p{z) for convenience. We will assume 

that the density and thus the potential are real, which 
effectively corresponds to a phase quenched approxima¬ 
tion. This approximation while limited in bulk 
still contains useful physical information at the edge of 
the Dirac spectrum [23] . 

Following [22], we substitute the collective potential 
(p{z) to the collective density p(z) through the Poisson 
equation 

VV(^) = -47r (^p{z) - (8) 

This change of variable involves the conformal anomaly 
in 2-dimensional curved space [52 [5S] 

Dp ^ det(-V2) Dp (9) 

with the Liouville action 

r2 = ^ / dzy/g + R{z)\n^) 

( 10 ) 

as we briefly detail in the Appendix. Here R{z) is the 
scalar Ricci curvature on the curved 2-dimensional com¬ 
plex manifold of eigenvalues. The determinant in ([^ 
induces a gravitational anomaly [22] . The final result for 
the effective action is 

Z[/3,p,iVy,m/]« (det(-V2))^ 

with the effective quenched action F = Do-l-ro-l-Fi-|-r 2 
and 

bo = ^ J dz^{{Vpf - R{z)p-AvB{z)p) 

Fi = ^ 0 J dzy/gp\np 

^^0 = ^ / dzdz'y/gg'hi\z- z'\‘^ (^Po{z) - 

( 12 ) 
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Here B{z) = B — TraS{z) and Fi = 0 for Dirac quarks 
with 1 / = 1/2. ( n][T2 ) differs from the one in [22] in three 
ways: 1/ Both p[z) and (p{z) are z-even because of chiral 
symmetry; 2/ i? is of order N-, 3/ B{z) carries a magnetic 
vortex which will be exploited below. We note that F is 
real and of order N'^. 

Although ( fnp^ was derived using the 2 -matrix 
model, we observe that each of its contributions are 
generic and therefore not specific to the 2 -matrix model. 
The exception is w{z) which is model specific and in this 
case quasi-harmonic. This suggests that (11) is the effec¬ 
tive partition function for QCD Dirac spectra at finite p, 
in the spontaneously broken phase provided that w{z) in 
Q is extended to include non-quasi-harmonic potentials. 


We note the quadratic form of the magnetic-like contri¬ 
bution, and the mixed curvature-magnetic contribution 
which is a Casimir effect. Recently, a similar mixed term 
between the gauge holonomy and the curvature in hyper¬ 
bolic space was noted in the cosmological context [25] . 
The quenched quark condensate in the ZMZ is (m/ = 

0 ) 

^ {m)o ~ Vg (*27rV^-5(z, z'))z'=o(16) 

with \7^Q{z,z') = S{z — z')f-y/g in curved eigenvalue 
space. In zeroth order in the curvature 2ttQ{z,z') = 
ln|z — z'\ and ([T^ vanishes 


IV. GROUND STATE 


In this section we analyze the ground state properties 
following from ([^ in terms of the collective potential in 
(11). We will explicit the solution to 5 F/i5:/j = 0 in the 
linearized approximation. The result is an expression for 
the mean charged density in the droplet po{z) without 
the contribution from the Vy-exponent in (11). We will 


refer to this solution as the quenched saddle point which 
is not to be confused with the standard quenched ap¬ 
proximation using in the litterature. We will use this 
result to analyze the mean contribution of the low-lying 
quark eigenmodes in the ZMZ to both the bulk energy 
and quark condensate at finite chemical potential and for 
fixed curvature i? « I/R 4 . Again, because of the generic 
nature of 0 as we noted earlier, we expect the results 
to reflect on the QCD eigenvalue spectrum at finite p 
beyond the 2 -matrix model. 


(17) 

as the last integral is zero when carried over the elliptic 
droplet ([^, 

L ^ “ l-l (^) = 0 

1 /2 

with a_(a;) = a_ (l — x^/a^) . This result is conform 

with the phase quenched limit, with the inside of the 
droplet breaking conformal symmetry with a condensa¬ 
tion of mixed pairs made of a quark and a conjugate 
quark [HUH]- 


B. Sub-leading 


A. Leading 

With this in mind, the quenched saddle point equation 
i5F /Sip = 0 in the linearized approximation yields the 
mean density 


, . I'B R va - . . va . 


(13) 


For vanishingly small p^, (13) yields 7 rpo(z) « N/Ap^ in 
agreement with the quenched asymptotic density in the 
weak non-hermiticity limit [24j . The quenched energy in 
the ZMZ is 


80 « -lnZ[/3,^,0,0] « J^Vg Po{z)W (z) (14) 

with pertinent changes in W(z) in curved space. Using 
([T^ we obtain 


» -T I „B 


+ )^^) /( VaM-) ( 15 ) 


The unquenched but subleading results follow from 
the linearized saddle point in (11) including the Nf- 
contribution. Specifically 


sr 


Sip{z 




The N'j contribution follows from the arguments pre¬ 
sented in [ 22 ] for the emergence of a conformal dimen¬ 
sion associated to in (11). Solving for (p in (18) 


yields the unquenched density to order 0{NfN°) 


vN f 

Pf{z) « Po{z) - — (S(z -f iruf) -k S(z - imf)) (19) 

The unquenched energy in the ZMZ in flat eigenvalue 
space is 


£■/ « -i ^ Pf{z) {W{z) - Nf\n{z^ + m})) (20) 


which differs from (14) to order 0{NfN) by 
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£f — So 


vNf I B 


ln(z^ 


■m)) 




( 21 ) 

The unquenched chiral condensate to the same order fol¬ 
lows from V 4 {qq)f = dSf/druf. We note that | (qq) f \ < 

I {qq)o |. Indeed, (21) yields to order NjN 


14 (( 99 )/ - ( 99 ) ( 


AT 

Nfrrif — 


a+TT 

a_|_ -I- a_ 


-f (1 t) 022) 


with uip = B/N and the rescaled momentum k = k/'/N. 
For /3 = 2, the screening length follows from kg + 24^;^ = 
0 or Zs « 1/^/N. Back to z-space 


5(5) 


(Pk 


r ~ - 

o-^i^-^Sik) = -^G 

TT 




(25) 


G(r) is Meijer G-function shown in Fig. with a loga¬ 
rithmic core G(r <?; Is) ~ —Inr, a hole of range Is and 
an asymptotic tail G(r Is) « . 


The effect of the fermion determinant on the 2- 
dimensional droplet amounts to inserting two static 
charges of vNf /2 at the mirror locations 2 = ztiirif. The 
charges are strongly screened as we now detail. 


V. SCREENING AND PLASMONS 

In this section we will explicit the screening nature of 
the droplet of eigenvalues viewed as a 2-dimensional 1- 
specie plasma with unit charges e = 1. This is one of 
the universal feature of the effective action 0 which 
goes beyond the 2-matrix model used. In particular, we 
will derive the static structure factor of the droplet us¬ 
ing a small and longitudinal deformation of the droplet 
density. We will observe that the static pole structure 
emerging from the structure factor nicely extends time¬ 
like to the plasmon pole contribution we have derived 
recently using a hydrodynamical analysis m. with an 
improved estimate for a relaxation time in QCD at finite 

fi. 


A. Static Structure Factor 



FIG. 2: Radial structure factor G(r) versus r. 


By rescaling z = ^/Nz in the microscopic limit, both 
B and po are of order 1 and the droplet fills out the entire 
2-plane in flat space. Within the center of the droplet, 
the (p(±zm^)-insertions do not affect the density and its 
response. The electro-static properties of the droplet are 
captured by the structure factor 


S{k) = 



(23) 


The longitudinal deformation 5p{z) « —poS/ ■ 4> yields 


S{k) « Pq/c^ ^|(/)(fc)py. It follows by expanding F in 
to quadratic order in 5p. The result is 


§ 


B. Relaxation Time 


The pole in (24) is the static limit of the longitudi¬ 
nal plasmon mode in the droplet. Indeed, ( [24| ) together 
with the hydrodynamical plasmon analysis in [19] implies 
the non-linear plasmon dispersion relation for the time- 
dependent longitudinal modes ipi^) —>■ 4 >{t, z) 


(y* + ^ $it, z) « 0 ( 26 ) 

which confirms and extends the hydrodynamical re¬ 
sult [12] 


Sik) 


UJp 


k‘^uj 


2 

p 




Aht 

48 ojp 




(24) 


'p 


4 


0(t,z)«O (27) 
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Now consider a time-dependent deviation of the eigen¬ 
value density through 5p{t, z) « —poV • 4>(t, z) in the 
complex 2-plane. For large droplets, the hydrodynam- 
ical arguments presented in [TS] show that the defor¬ 
mation relaxes through Euler equations with a large 


time asymptotics controlled by the plasmon branch (27), 
i.e. 5p[t —>■ oo,z)/po ~ with a relaxation time 

Tr ~ l/2ujp. If the initial condition for the eigenvalue 
distribution p(0, z) is chosen to describe a chirally sym¬ 
metric phase at finite p, then the time it takes for the 
distribution to relax to the spontaneously broken phase 
is again given by the plasmon branch. From (26) it fol¬ 


lows that the relaxation time for spontaneously break¬ 
ing/restoring chiral symmetry in QCD at finite chemical 
potential in droplets of finite sizes A = is now of 
order 


Tn 


1 

3-2 1 I /3 1 

4 .4 48 u>pA'^ 


(28) 


with « —1/A and Wp = B/N = 27^/3/A. Although 
( 26|28 ) follow from the model with /3 = 2, they nicely 
agree and extend the results in m following from the 
hydrodynamical arguments for the three Dyson ensem¬ 
bles /3 = 1,2,4. 


VI. BERRY PHASE 

In this section we will show that if the Dirac eigen¬ 
values where to change adiabatically with some mathe¬ 
matical time, i.e. Zi —>■ Zi(t), then the effective action 
will develop among other thinghs a geometrical contri¬ 
bution of the Wess-Zumino-Witten type. The origin of 
this term will be traced to a particular contribution in 
the eigenvalue measure which distinguishes between real 
(/3 = 1) or complex (/? = 2) representations for the 
quarks, i.e. whether the underlying quarks are Majorana 
or Dirac particles. This suggests that an adiabatically 
deformed droplet of eigenvalues at finite pL maps onto a 
2-dimensional fermionic system in a magnetic held. 

With this in mind, the magnetic held induced by the 
2-matrix model is B{z) = B — Tvad^z). The hrst con¬ 
tribution dehnes the mean density vB/2t: in hat space. 
The second contribution is a magnetic vortex of strength 
a/2 (z —>■ z) 


^ _ a £ij Zj 

O I 12 

z z 


(29) 


(29) is multivalued and generates a Berry phase 


For that, consider an adiabatic time-dependent change in 
the eigenvalue through z —> z{t) and v{t) = z{t)/\z{t)\. 
In the presence of (29), an anomalous Berry contribution 


is generated for each eigenvalue change as 


The form notation is subsumed. The line integral is over 
the mathematical time which is understood as a relax¬ 
ation time for Dirac spectra in l-|-4-dimensions [T5]. Each 
of the eigenvalue when adiabatically rotated around the 
origin accumulates a phase Si = aTr. For Dirac quarks 
in the complex representation and particle-anti-particle 
symmetry, /3 = 2 and a = 1. The phase accumulation is 

TT. 


To generalize (30) to all particles in the droplet we 
borrow from the argu ments in [30], and covariantize and 
localize u®(t) in (30) by using the embedding v^{t) —>■ 
v^{t, z) = (1, v{t, zB in 1 -I- 2 dimensions. Thus 


1 


1 


Sb = -olpqA / vdv —>■ -apii 


vdv 


(31) 


' 1 - 1-2 


which is a Wess-Zumino-Witten type term. 


VII. FERMI SURFACE 

The complex eigenvalues form a droplet in the com¬ 
plex 2-plane for all three Dyson ensembles. However, for 
13 = 2 with the underlying quarks in the Dirac representa¬ 
tion, adiabatically changing eigenvalues behave as mov¬ 
ing fermions in 2-dimensions with a fixed magnetic field 
B as defined in ([^ with filling fraction v = 1/(3 = 1/2 
since the mean droplet density is po = vB/2tt. Note 
that for (3 = 1 with the underlying quarks in the real 
or Majorana representation, the Berry phase is zero with 
no identification to fermions. We now explore the conse¬ 
quences of this relationship to fermions in 2-dimensions. 


A. Anomalous Fermi Surface 


We identify the Fermi momentum of the droplet as 
Po fc|^/47r or kp = \/2vB for (3 = 2. The anomalous 
transport on the Fermi surface follows by setting = 
kpv^ —>■ kF{vF,v(t, z)) in (31) and gauging by minimal 
substitution. The Fermi velocity vf < 1- Thus 


Sb = ^ / ^g{p +eA +sw)d{p +eA +sw) (32) 
87r JiF2 

Here A is a U(l) gauge held with F = dA and oj a U(l) 
spin connection with R = duj. Similar anomalous effec¬ 
tive actions for the quantum Hall effects were recently 
discussed in [221EJ. 


The U(l) current J = SSb/5A stemming from (32) is 
anomalous 




■/~9 


kpv^ 


8 tt 


F 


Stt 


R 


(33) 


Si= A* 


zdt=- 


zdz a f , 

2 , <“> 


The kp = 'j2vB contribution in (33) is the analogue of 
the chiral vortical effect in Fermi surfaces threaded by a 
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Berry phase in higher dimensions m- The anomalous 
contribution to the density is 


1 6 Sb 
= -tttt 

e dAQ[z) 


= 

47r 47r 


(34) 


A comparison of (34) with po in (13) suggests that e = 
2/a/3 and s = l/2a. For /3 = 2 and a = 1, we have 
e = 1 and s = 1/2. The droplet of eigenvalues maps onto 
a fermionic droplets in 2-dimensions with constituents of 
charge e = 1 and spin s = 1/2. 


mixed Chern-Simons term in the Fermi liqnid droplet is 
constructed following the arguments in |30j to which we 
refer to for more details. Since ea = 2//3 = 2^, the result 
in our case is 


5'n 


V 

47r 



d(se) 

e 


ef(e,s) 


x2 



Aosfl-h (A X V6» 


(36) 


B. Relation to Composite Fermions 


The Fermi distribution f(e, s) = l/(l-l-e('^ s/i_F)/r) 
composite fermions with particle-hole symmetry satisfies 


Recently, Son has argued that in a non-zero magnetic 
field the composite fermions of the fractional quantum 
Hall effect at half filling, have a finite density and live in 
a zero magnetic field |^. He argued that they exhibit 
particle-hole symmetry and that they are Dirac particles. 
Their ground state is a Fermi liquid with a Berry phase 
of TT. 

Our construction suggests that the anomalous trans¬ 
port of composite fermions in the fractional quantum Hall 
effect at half filling in 1 -|- 2 dimensions with a magnetic 
field, share similarities with a 2-matrix model of the ZMZ 
zone of Dirac quarks with /3 = 2 in l-|-4 dimensions in the 
spontaneously broken phase at finite chemical potential 
/i. The suggestion that composite fermions at half filling 
behave as Dirac particles with a Wess-Zumino-Witten 
term of the type (32) is falsifiable as it leads to specific 
and measurable predictions as we now detail. 


/ d{se) f (e, s) j = fj,p = kp = 'j2vB (37) 




',s=± ' 


X] / die){{e,s)\ =T\ii(2 + 2chi^ 


(36) yields an anomalous and D-driven contribution to 
the composite fermion density 


essentially from the Coriolis force contribution. The lim¬ 
iting results are measurable shifts in the density 


C. Rotating Fermi Surface 

We recall that novel chiral vortical effects were noted 
in rotating Weyl Fermi liquids with a Berry phase in 
1-1-3 dimensions [30] . We now show that similar effects 
are expected for composite fermions at filling fraction 
i/ = l/2 in 1-1-2 dimensions if they are Dirac particles. 
Throughout this sub-section B is understood as a real 
(residual) magnetic field. 

Consider now that the 2-dimensional fermions are 
in a rotating frame along the z-direction with veloc¬ 
ity n = riz. Each fermion experiences centrifugation 
that is best captured by the gravito-electro-magnetic 
fields Eg = p^VO and Bg = with a metric goo = 
1 — -I- 29 [32] • S acts as a gravitational-like po¬ 

tential. In the Fermi sea, the inertial 2-force on a quasi¬ 
particle of 3-momentum —>■ (e, skp) is the Lorentz-like 
force [30] 


Eg = eV6» + kpxsn (35) 

with the Coriolis force manifest. Here s = ± for a particle 
or anti-particle (hole). The corresponding Berry induced 


r ^ '^T\n 2 ^ 

Pfi\Z) —>■ -D 


T ^ kp = yj2i'B 


Pn(z) ^ T ^kp = V2i/B (39) 

47r 


The 0-contribution in (36) can be exploited using an 


observation by Luttinger who noted that the effect of a 
temperature gradient can be balanced by a gravitational 
potential |3Sj. The response of the composite fermions 
in a Fermi disc to a small temperature gradient can be 
captured by V0 = — VlnT. As a result, an anomalous 
and measurable stationary gradient flow develops 


dQ,i {z) 


SSn 

SA, 


u\/2vB 
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eijVjlnT 


(40) 


The results of this sub-section rely only on the interpre¬ 
tation that composite fermions are Dirac particles with 
particle-antiparticle (hole) symmetry and the origin of 
the Wess-Zumino-Witten term ( |3^ , with B a real (resid¬ 
ual) magnetic field as we noted earlier. Therefore they 
are more general than the correspondence we have devel¬ 
oped. 

















VIII. CONCLUSIONS 

The QCD Dirac spectrum at finite chemical potential 
following from a 2-matrix model at large N maps onto 
a generic but anomalous effective theory in 2-dimensions 
with both gauge and conformal anomalies in bulk, much 
like in the fractional quantum Hall effect [22]. The com¬ 
plex eigenvalues form a droplet threaded by a magnetic 
vortex at the center. Adiabatic changes in the Dirac 
eigenvalues are characterized by a Berry phase that leads 
to a Wess-Zumino-Witten term. The anomalous trans¬ 
port of composite fermions at half hlling in 1 -|- 2 dimen¬ 
sions share similarities with a 2-matrix model of Dirac 
particles in 1 -|- 4 dimensions in the spontaneously broken 
phase at finite chemical potential. 

IX. ACKNOWLEDGEMENTS 

We thank Sasha Abanov and Piotr Warchol for dis¬ 
cussions. This work is supported in part by the U.S. 
Department of Energy under Contracts No. DE-FG- 
88ER40388. 


X. APPENDIX: CONFORMAL ANOMALY 


Appendix, we give a brief account of the derivation fol¬ 
lowing [22]. The change from the collective density p to 
the collective potential (p follows formally from Q as 


deti-V"^) Dip = ( / Dp (41) 


with ? 7 , fj as Grassmanians. The regularized measure 
DrjDfjDp depends implicitly on the collective density p 
and thus the collective potential p through regulariza¬ 
tion [22l|28]. As a result, p can be treated as a conformal 
factor of the mathematical metric on the eigenvalues, 
and removed from the measure by a conformal transfor¬ 
mation of the coordinates 


dV 2 = y/gdz ^ dV 2 /p (42) 


The emergence of the Liouville effective action in (|^ 
10) follows from the conformal anomaly [221[2S1- In this 


The result is (9][l0) after using the central charges c^p = 


-1-1 and Cri = Cfj = —1. 
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